For submission to MNRAS. Not for distribution 



The Formation and Evolution of Virgo Cluster Galaxies - 
I. Broadband Optical & Infrared Colours 

Joel C. Roediger and Stephane Courteau 

Department of Physics, Engineering Physics & Astronomy, Queen 's University, Kingston, 

Ontario, Canada 

Michael McDonald 

Dept. of Astronomy, University of Maryland, College Park, MD 

and 

Lauren A. MacArthur 

Herzberg Institute of Astrophysics, National Research Council of Canada, 5071 West Saanich 

Road, Victoria, BC, Canada 

j roediger, courteau@astro. queen su . ca, mcdonald@ astro . umd . edu, 
lauren . macarthur@nrc-cnrc . gc . ca 

ABSTRACT 

We use a combination of deep optical (gri) and near-infrared (H) photometry to 
study the radially-resolved colours of a broad sample of Virgo cluster galaxies. For 
most galaxy types, we find that the median g-H colour gradient is either flat (gas-poor 
giants and gas-rich dwarfs) or negative (i.e., colours become bluer with increasing 
radius; gas-poor dwarfs, spirals, and gas-poor peculiars). Also, later-type galaxies 
typically exhibit more negative gradients than early-types. Given the lack of a correla- 
tion between the central colours and axis ratios of Virgo spiral galaxies, we argue that 
dust likely plays a small role, if at all, in setting those colour gradients. We search for 
possible correlations between galaxy colour and photometric structure or environment 
and find that the Virgo galaxy colours become redder with increasing concentration, 
luminosity and surface brightness, while no dependence with cluster-centric radius 
or local galaxy density is detected (over a range of ~2 Mpc and ~ 3 — 16 Mpc-2, 
respectively). However, the colours of gas-rich Virgo galaxies do correlate with neu- 
tral gas deficiency, such that these galaxies become redder with higher deficiencies. 
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Comparisons with stellar population models suggest that these colour gradients arise 
principally from variations in stellar metallicity within these galaxies, while age vari- 
ations only make a significant contribution to the colour gradients of Virgo irregulars. 
This paper serves as an introduction to a detailed stellar population analysis of Virgo 
cluster galaxies by the same authors (Roediger et al 201 lb). 

Subject headings: galaxies: clusters: individual (Virgo) — galaxies: dwarf — galax- 
ies: elliptical — galaxies: lenticular — galaxies: evolution — galaxies: formation — 
galaxies: irregular — galaxies: peculiar — galaxies: spiral — galaxies: stellar content 
— galaxies: structure 



1. INTRODUCTION 



The broadband colours of galaxies reveal information about the ages and metallicities of their 
stellar populations (as well as their dust content). Galaxy colours thus provide a fundamental 
means for the investigation of both their star formation h istories and chem i cal evolution, a fact 



which is reflected by the rich history of this field (e.g.. IVader et al.l 11 988 



Goudfro oii et al]ll994t Ide Jondll996l: Bell & de Jon g 2000; 



2004 iTavloret all 120051) " 

since the pioneering works of 

d 1976b . 



Peletier et al. 1990; 



MacArthur et al 



Searle et all (119731) and 



2004 



Ivan Zee et al. 



Tinslev & Gunn 



The alternative to colour-based analyses of galaxies' stellar populations involv es a spectro 



scopic approach, either through the measurement of absorption line strengths (e.g., Faber 



Terlevich et al 



2000; 



198lllBurstein et al.ll984llwbrfhey et al.l 1 9 921: llCobavashi & Arimotol 1999c 
Gallazzi et al.l2005l : |Peletier et al.ll2007|) or the synthesis of entire spectra (e.g., 



1973 



Trageretal 



MacArthur et al. 



2009). The benefits of the s pectroscopic a pproach include a cleaner separation of the correlated 
age and metallicity effects (Worthev 1994) and, in the case of the Lick indices, a global insensi- 
tivity to dust reddening (|MacArthurll2005|) . Spectroscopic stellar population studies demand some 
appreciable trade-offs compared to colour-based analyses however, such as smaller sample sizes 
and shallower depths. Colours then offer the most efficient means to simultaneously obtain deep 
radially-resolved and representative stellar population information for large numbers of galaxies. 

Through comparison with stellar population models, galaxies' colour gradients may be de- 
composed into corresponding gradients of age and metallicity. The latter data provide a record of 
the build-up of galaxies' stellar contents and so are naturally of great value to galaxy formation and 
evolution models. In modelli ng galaxy colours, a wide and well-sampled wavelength baseline (i.e., 
ultraviolet, UV, to infrared; iGavazzi et al.ll2002|) is certainly desirable, not only to better constrain 
the fits, but also to overcome the age-metallicity degeneracy which plagues optical stellar popula- 
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tion studies ( Wortheylll994|) . Although the availability of space-based imagin g (i.e., GALEX) ha s 
provided a much desired wavelength expansion of colour-based studies (e.g.. lBoselli et al.l l2008). 
the use of optical and ne ar-infrared (NIR) data alone already goes a long way towards overcoming 
the above degeneracy (|Cardiel et al.ll2003|) . 



ies (E/SO), the consensus seems to be that they possess negative colour gradients 


(Peletier & Valentiin 


1989; 


Peletier et al 


. 1990: 


Goudfrooii et al. 


1994; 


Silva & Elston 


1994 


. Michard 1999: Wu et al. 


2005; 


James et al.l2006c La Barbera & de Carvalho 


2009 


; Gonzalez-Perez et al. 


2011 


, La Barbera et al. 


2010), the cause of which is popularly attributed to an outward decrease in stellar metallicity 


(Vaderetal.lll988: 


Franx et al.lll989tlTamura et al.lboOO: 


Tamura & OhtalboOO. 


2003; 


La Barbera et al. 


2004) 


. Positive UV and optical colour gradients found in both local and high-z E/SO's dMichardl 


2000; 


Moth & Elstonbooi 


Menanteau et al 


2004 


Ferreras et al. 


2009; 


Suh et al.ll2010). however. 



imply that centralized star formation has recently occurred, thereby arguing for at least a partial 
age gradient contribution. 

Discerning the nature of spiral galaxies' colour gradients is in herently complicated due to the 
presence of dust, gas, and multiple stellar populations within them ( Kim & An"nlll99C ; Gadotti & dos Anjos 
200iI : Ikoo et alJbooiFfavlor et all2005l : lGanda et all2009l : lLiu et al.ll2Q09h . This complexity may 
be lifted though by distinguishing between early- and late-type spirals , wherein the former have 
flat or negative colou r gradients (ITerndrup et al.l 1 19941: Ide Jongl 1 19961 : IPeletier & Balcellslll996 : 
Moriondo et al. 2001 ; Ganda et al. 2009 ) while the latter are described by positive colour gra- 
dients, but with a wide dispersion (|Kim & Annlll990l : iTaylor et al.l 120051 : iTamm & Tenjesll2006l : 
Liu et al.ll2009|) . Given such a variety, it is not surprising that spiral galaxies' colour gradients 
have often been ascribed to both age and metallicity effects, with dust making a second-order or 



20011 : iMacArthur et al.ll2004l ; however, see 



altogether negligible contrib ution (|de Jon g 1996; Peletier & B alcellsl 1 1 996c iGadotti & dos Anjos 



Ganda et al. 



2009 



Lastly, the colour gradients amongst the remaining galaxy types are more uniform than those 
found in spirals. That is, both peculiar and (gas-rich and gas-poor) dwarf galaxies are typically 



James 



described by positive colour gradients (IVader et . al.lll988t Uerieneta 



1994; Pa paderos et all 1 19961 : iTullv et all 1 19961 : Barazza et all 



2000; Kim & Ann 



2003; 



van Zee et al. 



1990; 



2004; 



Taylor et al.ll2005|) . Since the notion of such galaxies having a higher metallicity in their outskirts 
(where the potential is weaker) seems highly counter-intuitive, these colour g radients have com - 
monly been attributed to underlying positive age gradients within therrQ (e.g., Vader et al. 1988). 
Negative metallicity gradients, on the other hand, may very well contribute to the flat and negative 



The effects of dust and recent star formation may, however, also contribute to the positive colour gradients in 
gas-poor dwarfs. 
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colour gradients that have been observed within a minority of dw arf galaxies ([Patterson & Thuan 
199d : IParodi et allbooi barazza et all2003l: Ivan Zee et al1b004h . 



Aside from the typical colour gradient for each principal morphological type, the structural or 
environmental parameters which correlate with galaxies' colours and colour gradients are also of 
interest to understanding galaxy for mation. The most prevalent of such correlations for all galaxy 
types are those against luminosity (IBalcells & Peletier 1994 : Ijansen et al.l 12000: Tamura & Ohta 



2003 



Munoz-Mateos et al. 



2007; 



Mac Arthur et alJ2004HWest et al. 



1996;Pierini 



Liu et al.1 



2009 



Tortora et a 



2002; Gadot ti & dos Anios 200 1 ; 



1999; Parodi et al. 2()oJ;|Tamura & Ohta 2003 



Ferreras et al 



20 lX mass dVaderetal. 



2005 



1988 



2009l:lLa Barbera et al.lboioLISuh et alhoiol) , col our dTullv et a 



2009k . and size (bereyall et al 



Munoz-Mat eos et al.ll2007l : iLiu et al. 



2009: lTortora et al 



2010() . in the sense that more luminous, more massive, redder, and larger galaxies have reddest 
colours and flattest colour gradients. In additio n, spiral galaxies have bluer colours and steeper 
colou r gradients for low er surface brightnesses (|de Jonglll996l : ILiu et alj|2009l) and larger Hi con- 
tents (IWest et al.ll2009|) . Less uniformity has been found amongst galaxy types with respect to the 
environmental dependen ce of their co l our information, however. While the colour gradients of 
gas-poo r giants flatten dK o & Im 120051 : La Barbera et al. 2005) and the colours of gas-rich dwarfs 
redden (|Gallagher & Hunterl 11986c IParodi et al.l 120021) towards high densities, no correspo nding 
trends have been found amongst spiral galaxies ( Koo et al.1120051: [Munoz-Mateos et al.ll2007|) . 



Despite the seemingly coherent picture presented above based upon piecemeal information, 
a homogeneous analysis of galaxy colours and colour gradients across all galaxy types which 
incorporates NIR data would enable a valuable confirmation of that picture. Thus, we present here 
the resolved optical and optical-NIR colours of a sample of Virgo cluster galaxies that spans the 
full range in galaxy morphology. This paper is organized as follows. In Owe discuss the salient 
aspects and organisation of our photometric database. Broadband colour gradients, as well as 
morphological, structural, and environmental trends with galaxy colours are shown and discussed 
in the context of previous work in §(3] The interpretation of our results in terms of the stellar 
populations of Virgo cluster galaxies, however, is addressed in a companion paper (Roediger et al. 
2011b; hereafter, Paper II). Conclusions are provi ded in §|4l Thro ughout this paper, we assume a 
distance of 16.5 Mpc to all Virgo cluster galaxies (|Mei et al.ll2007|) and refer to both SO and SO/a 
galaxy types simply as "SO". 



2. DATA 

A systematic study of the colours and stellar populations within galaxies of all types demands 
as complete a volume-limited sample as possible. Given that the Virgo cluster is the closest large 
concentration of galaxies whose members span the full range of galaxian parameter space, we use 
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for our study the radially-resolved griH photometric database derived from the imaging survey of 
Virgo galaxies of McDonald et al. (201 1, hereafter Mil). Only those aspects of Mil which are 
relevant to our study are discussed below; further details can be found in Ml 1. 



The Mil sample includes 285 galaxies selected from the Virgo Cluster Catalog (|Binggeli et al. 

1985L hereafter VCC ); it is nearly complete down to B ~ 16 (M B ~ -15). For each galaxy, Sloan 



Digital Sky Survey (lAdelman-McCarthv et all I2007L I2008L hereafter SDSS) gri- and separate 
//-band imaging is available. The spatial distribution and completeness of this sample are shown 
in Figures 1 & 2 o f Mil. The Mil sample extends out to the 6° radius of the Virgo cluster 
(|McLaughlinl 1 1999Q . such that the sample spans a wide range in environmental conditions (i.e., 
galaxy densities S ~ 3-16 Mpc" 2 and gas deficiencies Defm ~ -0.4-2.3 dex; see 32.31 As seen 
in Mil, a natural bias for gas-poor galaxies exists in the sample given the morphology-density 
relation (lDresslerlll980|) . Nevertheless, the Mil sample contains 86 gas-rich galaxies (56 giants 
and 30 dwarfs; see Tabled), thereby allowing for meaningful comparisons of the colours, and thus 
stellar populations, between major Virgo galaxy types. 



2.1. Photometry 



Radially-resolved optical (ugriz) photometry for the Mil galaxies was extracted from SDSS 
imaging. However, the u- and z-bands were excluded from our analysis in light of their low signal- 
to-noise (hereafter S/N) and the significant number of light profiles having obvious sky subtraction 
errors at those wavelengths (Mil). The Mil NIR photo metry came from a mixture of imaging 
sources, including the GOLD Mine (iGavazzi et al. 2003 , 78 galaxies) and Two Micron All-Sky 
Survey (ISkrutskie et al.l |2006[ 20 galaxies) databases, while //-band imaging for the remaining 
(187) galaxies was obtained by Mil from April 2005 to May 2008 on Mauna Kea telescopes 
[ULBCAM (UH 2.2 m); WIRCAM (CFHT); WFCAM (UKIRT)]. The spatial extent of the NIR 
light profiles is, on average, ~60% that of the corresponding /-band profiles for ~200 galaxies (see 
§12.41 for our definition of profile extent), while the optical and NIR profiles reach median surface 
brightness levels of r ~ 26.5 mag arcsec" 2 and H ~ 23 mag arcsec" 2 , respectively. The high quality 
of the Mil light profiles enables a reliable investigation of Virgo galaxies' colours out to typically 
large galactocentric radii. 
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2.2. Surface brightness profiles and structural diagnostics 

The griH major-axis surface brightness (hereafter SB) profiles of Virgo galaxies were ex- 
tracted by Ml 1 with the astronomical image-processing software XVISTa|. To ensure that colour 
gradients are computed homogeneously, the SDSS SB profiles were degraded to the (poorer) reso- 
lution of the NIR profiles. 

We hav e also drawn upon the collection of non-parametric, //-band structural parameters 



presented in Mc Donald et al.l (|2009() to search for possible correlations between galaxy colour and 
structure. The structural parameters of interest include effective surface brightnesses (fj, e ), effective 
radii (r e ), concentrations (C 2 g) and apparent magnitudes (m H ). These non-parametric diagnostics 
all involve the extrapolation and integration of galaxies' NIR light profiles to infinity. We use the 
NIR values of the above structural diagnostics for their lower sensitivity to dust attenuation and 
greater sensitivity to the galaxy stellar masses. 



2.3. Environmental diagnostics 



To probe for an environmental dependence in the Virgo galaxy colours, we characterize the 
impact of environment with the neutral hydrogen deficiency parameter, Def HI . For a cluster galaxy 
with a given gas mass M HI , Def HI = (\ogM HI (D optJ )) - \ogM m , where (\ogM HI (D optJ )) is the typ- 
ical Hi mass for field galaxies of the same Hubble type, T, and optical diameter, D opt , such that 
positive Def H i values correspond to a dearth of neutral gas. The advantage of Def HI as a mea- 
sure of galaxy en vironment is its ind ependence of projection effects that plague other traditional 
(2D) diagnostics (|Solanes et al.1 l2002l hereafter S02). For completeness, we still appeal to the 
(projected) cluster-centric distances, D M %i, and local surface densities of the Mil galaxies, E, to 
assist our study of environmental effects on their colours. Cluster-centric distances were computed 



relativ e to M87 w hile assuming a common distance to all VCC galaxies of 16.5 Mpc (IMei et al. 



2007k We follow bresslerl (1980) in computing £ as the number density comprised by the ten 
nearest VCC neighbours to a given galaxy's position. 



We obtained De fni measurements fo r Mil galaxies from S02, iKoopmann & Kenneyl ((2004, 
hereafter KK04), and lGavazzi et alj (120051 . hereafter G05). The latter compilation is the most com- 
plete, providing Hi mass measurements for 91 of the 115 gas-rich galaxies in the Mil sample, as 
well as an additional 10 E/S0 systems. The other Hi mass compilations have Def m determinations 
for 28 (KK04) and 53 (S02) galaxies, only four of which are not included in G05. We compare 
Defm measurements in Figure [1] for galaxies in common for KK04 and G02 with S02, with the 



2 http://ganymede.nmsu.edu/holtz/xvista/ 
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red line in each plot denoting equality. While there is good consistency between the KK04 and 
S02 databases (r = 0.97; a = 0.11 dex), the measurements of G05 are systematically higher on 
average than the others, regardless of morphology. Presuming that the G05 measurements are at 
least internally consistent, we rely strictly on that database for Def HI values. 



2.4. Surface brightness profile preparations 

Prior to measuring Virgo galaxy colours and stellar populations, the Ml 1 SB profiles require 
special attention, which we now describe. The computation of accurate galaxy colour profiles 
relies on accounting for pixels that are spoiled by either foreground stars, scattered light, or low 
S/N. To correct for this, we first truncated the Mil SB profiles by searching outward in each of 
them for a sequence of adjacent isophotes with SB errors in excess of a given threshold a^ t h- The 
isophotes beyond that point were flagged as unreliable and their location marked as the truncation 
radius. 

In the truncation process, we used error thresholds of 0.25 and 0.20 mag arcsec" 2 for the Mil 
optical and NIR SB profiles, respectively. These values correspond to locations in the SB profiles 
where the point-to-point scatter becomes comparable to the statistical errors per isophote. These 
thresholds are also roughly equal to 0.1-0.2% of the sky level (or less) in each profile. 

We have also applied an adaptive radial binning scheme to the Mil S B profiles to increase the 



S/N per radial bin. While previous studies of galaxy colour gradients (e.g. jMacArthur et al.ll2004j) 
used an incremental radial binning scheme, ours requires that a minimum signal-to-noise S/N m m 
be met in all four bands in a given bin. Since an incremental scheme makes no such stipulation and 
typically has larger radial bins, the colour profiles of individual galaxies will then tend to extend 
further in this scheme than those determined from our adaptive scheme. However, these increased 
depths are overshadowed by the higher colour dispersions per radial bin and the inclusion of low 
S/N data points in galaxies' outskirts, thereby decreasing the quality of the colour profiles overall. 

Accounting only for the Poisson noise of the source and sky, we calculate the S/N per isophote 
according to, 



S S. 



gal 



N y/Nla+NZy 

_ /^gal x Qll 

\/ /igal X Cell + A*sky x Cell 

^" /ck, (i) 



V /%al A^sky 
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where yU ga i and fi sky are the isophote's SB and sky background, respectively, and C e u is its circum- 
ference. 

The radial binning scheme for each galaxy is based on the SB profile with the lowest overall 
S/N; the latter usually corresponds to the g-band, followed by the //-band. We enforce a minimum 
bin size for each galaxy, which is set by the maximum seeing disk amongst all bands and most 
often corresponds to the //-band value. The SB of a bin is determined by a weighted average of 
the isophotal SBs that it contains, while the bin location is set to its midpoint. Through extensive 
testing, we have found that using the cr^-fi median trends determined from the original photometry 
provided the most reliable way of estimating the SB error per bin. Lastly, we have carried out our 
analysis for the cases of S/N m i n = 10, 20, 30, 40 and 50 per bin. 



2.5. Sky errors 



The determination of accurate colour gradients for a galaxy demands reliable subtraction of 
the sky's contribution to the flux measured within each filter. An under- or over-estimated sky in a 
given band will artificially increase or decrease the corresponding SB measured in a galaxy's out- 
skirts, respectively, thereby skew ing its colour profile s. Sky errors dominate all other error sources 



in the study of galaxy colours (|Bell & de Jong 



2000). While Mil took great care in assessing sky 



level corrections, our profile truncation algorithm (§ 12.4) is already effective before the sky error 
envelopes become prohibitively large (i.e., cr sky ~ 0.2% of the sky level). We revisit the issue of 
sky errors in § 13.11 



3. RESULTS & DISCUSSION 

We now consider the broadband colour profiles of Virgo galaxies, as well as seek out possible 
morphological, structural and environmental trends in their central colours. These data offer a 
view of stellar populations both within and amongst galaxies, free of any assumptions about their 
underlying star formation histories and chemical evolution. 



3.1. Colour profiles and gradients 

The wide frequency baseline of the g-H colour provides us with our best constraint on the 
broad spectral shapes of Virgo galaxies, and thus on their stellar populations. The median g-H 
profiles and rms dispersions for Virgo galaxies of different morphologies are shown in Figure |2] 
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by the solid and dashed lines, respectively. For this exercise, the g-H colour profile for each 
galaxy was created by subtracting its binned //-band light profile from its binned g-band light 
profile. We also scaled each colour profile according to the NIR effective radius of the galaxy, 
thereby removing the size variation that exists amongst galaxies of various Hubble types. Given 
our imposed S/N constraints, we find that the median colour profiles for the many Virgo galaxy 
types span a large range, extending out to ~6 r e for the gas-poor giants, down to >2 r e for the 
gas-rich dwarfs and irregulars. 

As a complement to Figure |2l we list in Table \T\ the relevant colour information for each 
morphological group (column 1). Columns 2-4 contain the median g-H colour gradients, where 
the gradient for each galaxy was measured from a linear fit to the g-H profile. The gradients were 
computed for three choices of radial coordinate, so that the median gradients could be expressed 
either in terms of scaled (r e ; column 2) or physical radii (the latter were computed on both a linear 
and logarithmic scale; columns 3-4). Column 5 shows the median g-H central colours, while 
column 6 provides the number of galaxies in each morphological group. Both the central colours 
and colour gradients for each of our galaxies, along with their surface brightness profiles and bulge- 
disk decompositions, are available at the Centre de Donnees astronomiques de Strasbourg and our 
own website§. 

In fitting galaxies' colour profiles, we neglected the first radial bin in order to avoid adverse 
PSF blurring effects from skewing our computed gradients. We show in Figure |3] the size of the 
first radial bin versus the FWHM of the maximum seeing disk for all Virgo galaxies. The proximity 
of many galaxies to the line of equality (red) indicates that seeing effects within the first bin of the 
colour profiles for those systems are likely non-negligible. 

Other effects which could potentially skew the median colour gradient estimates provided in 
Table \T\ include the use of a variable radial range in individual colour profile fits (given our choice 
of binning scheme) and of different S/N min values in adaptively binning the Mil SB profiles. As 
shown in Figures @] and |5] though, neither of these two effects presents a significant risk towards 
the robustness of our colour gradient determinations. 

Examination of Table \T\ reveals a large dispersion in the colour gradient distribution of each 
Virgo galaxy type. These dispersions may be either intrinsic to Virgo cluster galaxies, or largely the 
result of measurement errors in their SB profiles. As mentioned in §12.5 L imperfect sky subtraction 
is the largest source of error in the computation of galaxy colour gradients. We can estimate the 
impact of sky errors on those gradients through the ± 1 a s ^ error envelopes for each M 1 1 SB profile, 
where a sky is the standard deviation in the sky levels measured in boxes along the periphery of a 
given Virgo galaxy image. We assess the colour gradient error by calculating the extrema of this 



www.astro.queensu.ca/virgo 
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quantity based on sky error envelopes. The distribution of errors in the g-H colour gradients of 
Virgo galaxies due to imperfect sky subtractions is shown in Figure [6] Comparing Figure [6] with 
Tabled] it can be seen that median sky errors contribute no more than ~16% of the dispersion 
found in the colour gradient distribution of each galaxy type; those dispersions are thus largely 
intrinsic to each galaxy population. 

It is reassuring that the distribution of the median colour gradients listed in Table \T\ accurately 
reflects our impressions from Figure |2] which we now describe^. Collectively, the median g-H 
profiles of Virgo galaxies exhibit intriguing features. First, most galaxy types (dE/dSO, Sa-Sd, 
Im, ?) are described by a negative colour gradient, albeit of differing strengths. The gas-poor 
E/SO's, Sdm+Sm spirals, star-forming BCDs, and morphologically-peculiar S? galaxy types, on 
the other hand, have either quasi-flat (E/SO, Sdm+Sm, BCD) or positive (S?) colour gradients. 
This variation in the g-H gradients of Virgo galaxies seems to correlate with morp hology, such 



that the gradient is more negative for later Hubble type s, as seen elsewhere (e.g., IJansen et al 



2000 ; La Barbera et al.ll2002 ; Gonzalez-Perez et ap201 1 ). Although the large dispersions reported 
in Table \T\ might draw suspicion on such a claim, comparisions of the gradient distributions for 
neighbouring morphologies via the Kolmogorov-Smirnov (KS) test confirms that the probability 
that any two distributions come from the same parent does not exceed ~40%. 

On its own, our detected morphological trend in the median colour gradients of Virgo galax- 
ies might seem surprising as previous studies of galaxy colours have failed to reveal the exis- 
tence of such gradient (i.e . , positive/negative) for most galaxy types (e . g.. Kim & Ann! Il99d ; 



Gadotti & dos Anjosl 1200 ll : iGrant et all 120051: Bunter & Elmegreen 



2006; Liu et al. 



2009). The 



dispersions in Table Q] show, however, that the colour gradient distribution for each Virgo galaxy 
type encompasses both positive and negative values, confirming that the variety in colour gradient 
sign found previou sly is reproduced here. Late-type spirals, on the other hand, have been shown by 
iTavlor et al.ll2Q05h as having a robust (positive) colour gradient sign, which, interestingly enough, 
is contrary to the typically negative gradient that we find for Virgo Scd-Sm spirals. This dis- 
crepancy may be attributed to the fact that Taylor et al. used UV-optical colours, which are quite 
sensitive to the ongoing star formation that is likely occurring within these gas-rich galaxies. 

Also interesting from Table [T]is the variation of the colour gradient dispersion amongst major 
galaxy types. For instance, the (rms) dispersion (in r e space) for the Virgo gas-poor giants and 
spirals is ~0.21 and ~0.36, respectively. The irregulars have intermediate values. This result 
agrees with the significa nt scatter in both the colours and colour gradients of spiral galaxies that has 
been reported by many (IPeletier & Balcellslll996t [Taylor et al .1120051 : iMunoz-Mateos et al] 120071 : 



The apparent contradictory behaviour exhibited in the outskirts of some median colour profiles (e.g., dE) is re- 
solved by considering that the colour data in the outskirts have lower S/N than that which comprises the interiors 
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Ganda et al J 12009c IWest et al.ll2009f) . The large colour ranges spanned by both Virgo spirals and 
irregulars may reflect a large variance in their stellar population properties due to their ongoing 
star formation. Dust attenuation may also be affecting these galaxies' colours, especially since the 
Mil sample contains a few unrealistically red (g-H ~ 4-5) Sa-Sd galaxies (see Figure |9). Whether 
the role of dust in setting galaxy colours is statistically important or not will be revisited in § 13.2.41 

In addition to the collective properties and trends noted above, inspection of the median colour 
profiles amongst major Virgo galaxy types reveals other noteworthy results. For instance, within 
their respective errors, the profiles for the two gas-poor giant types (E/SO) are remarkably similar; 
both exhibit a relatively constant colour (g-H ~ 3.2) out to ~4 r e , beyond which they become 
redder with radius. This suggests that Virgo E/SO's possess similar stellar populations and likely 
share a similar (extended) formation history; this idea is explored further in Paper II. On the other 
hand, Figure |2] and Table \T\ remind us that Virgo gas-poor dwarfs are systematically bluer and 
have more negative colour gradients than the gas-poor giants. Indeed, KS tests of the colour 
gradient distributions for Virgo gas-poor dwarfs and giants only have a ~7% probability of having a 
common parent distribution. The colour gradients that we have found for Virgo gas-po or dwarfs are 



consistent with these gala xies being less enriched overall than Virgo gas-poor giants (|Zaritsky et al. 



1994 iTremonti et al.ll2004c Paper II). 



3.1.1. Comparison with literature 



It is instru ctive to compare our opt ical-NIR colour gradient estimates with literature values for 
other galaxies. La Barbera et al.l (|2010|) measured d(^-//)/dlog(r/r £ .) colour gradients in 4546 giant 
gas-poor galaxies that overlap with both the SDSS and 2MASS to find a strongly negative mean 
value and large dispersion for them: -0.286 ± 0.293 mag/dex. The combined g-H colour gradient 
distribution for our Virgo E and SO galaxies, equal to -0.245 ±0.426 mag/dex, agrees well with 
La Barbera's result. Moreover, the analysi s of the g-K colo ur gradient distribution for a smaller 
sample of nearby giant gas-poor galaxies by IWu et al.l (|2005|) has also yielded a mean value similar 
to ours (-0.29 mag/dex). 



With respect to gas-rich dwarf galaxies, Hunter & Elmegreenl (|2006|) analysed the V-J colour 
gradients in both Magellanic irregulars and blue compact dwarfs and found that the mean colour 
gradient for these two galaxy types is -0.07 and +0.90 mag/kpc, respectively. While the mean g-H 
colour gradient that we have estimated for Virgo Im galaxy types (-0.1 1 mag/kpc) agrees well with 
the former, that for Virgo BCDs (+0. 1 1 mag/kpc) is much flatter. This discrepancy is reduced when 
we consider the fact that the standard deviations about both mean values is large (1.129 and 0.318 
mag/kpc, respectively) and that both distributions suffer from small sample sizes (N = 5 and 8, 
respectively). 
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Lastly, iMunoz-Mateos et al.l (I2007|) studied the FUV-K colour gradients for spiral galaxies 
of all types. Although a quantitative comparison of FUV-K and g-H colours and colour gradients 
seems ill-conceived for any galaxy type, it is interesting to note that the median FUV-K colour 
gradient in these authors' sample decreases, albeit weakly, from early- to late-type spirals (-0.015 
mag kpc" 1 for Sa-Sb's to -0.040 mag kpc" 1 for Scd+Sd's); their Sdm+Sm spirals do not follow 
this trend, however, having a nearly flat median gradient of -0.005 mag kpc" 1 . This morphological 
trend in the Munoz-Mateos et al. data is similar to the one which we have found in terms of the 
g-H colour gradients of Virgo cluster spirals. 

Based on these comparisons, encompassing a wide range of galaxy types, Virgo cluster galax- 
ies exhibit radial colour distributions that are apparently similar to those of field and cluster galaxies 
found in other nearby galaxy catalogs. 



3.2. Colour trends 



3.2.1. Morphology 



We now investigate trends in Virgo galaxy colours as a function of galaxy morphology, struc- 
ture, environment and axis ratio. Beginning with morphology, we show colour-colour diagrams in 
Figures [7]|9] for several galaxy types against the colour predictions from a stellar population model 
(black lines) in order to diagnose the effects of mean age (A) (g-r; dashed) and metallicity Z (r-H; 
solid) in setting the Virgo galaxy colours. The model is an exponential star formation history with 
a variable timescale r and several fixed metallicities. Nominal values of both the (A) and Z co- 
ordinates are indicated alongside the grid in each plot. Further details about this and other stellar 
population models are given in Paper II. The orange-red arrow represents the reddening vector for 
a dust screen, also discussed more fully in Paper II. Median values of both the calibration and sky 
errors in the Mil sample are shown in the lower-right comer of each plot, where the effect of the 
latter has been computed at both the centers and outskirts of Virgo galaxies. Note that both the 
model predictions and reddening vector discussed above have been included in Figures [10lfT2l as 
well. 

Figures [7]|9] show that the different Virgo galaxy types have unique colour profile combina- 
tions, and thus stellar contents. The colour profiles of the gas-poor systems (Figures [7]|8]) suggest 
that their stellar populations are old (>9 Gyr) and span a wi de range in metall i city; that is, these 



colour gradients are largely explain ed as a metallicity effect (|Vader et al.lll988E iFranx et al.lll989 



Tamura et al. 2000; Wu et al 



contribution (La Barbera et al 



2005), while age vari ations make a secondary, but n on-negligible, 



2010[ISuhetalJ20ld but see bothun & Greggl[l99ol) . The system 



atic differences in the median colour profiles between the Virgo gas-poor giants and dwarfs (§ |3.1t 
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are similarly explained by the stellar populations of the former being more metal-rich ( AlogZ < 
-0.25), but only slightly older (~0.5 Gyr), than the latter. Interestingly, the dE's and E's together 
define the metallicity extrema for the entire Virgo gas-poor galaxy population, as opposed to the 
dSO's and SO's. 

The Virgo irregulars (Figure |8) show the greatest colour range amongst all of the major galaxy 
types in this cluster such that, collectively, their stellar populations go from being old and metal- 
rich to young and metal-poor. This extreme colour variation extends to the Virgo irregulars as 
well, as their median colour profiles show that significant radial increases and decreases in age and 
metallicity, respectively, occur within each galaxy type. Such broad population gradients make it 
difficult to assign typical ages and metallicities to irregular galaxies. The offsetting effect that such 
radial behaviours in stellar populations have on colours helps explain previous reports of flat (or 



weak l y positive) optical colour gradients in irregular galaxies (IKim & Annlll99Q : |Papaderos et al. 



1996 



Taylor et al.l 120051). Moreover, th e redder disks that have been found in dlrr's relative to 



those in BCDs (|Patterson & Thuanl ll996) may be understood by the former having older outskirts. 



Age and metallicity gradients are also apparent within Virgo spirals (Figure |9), but with the 
latter being more significant, on average, than the former. Age and metallicity variations h ave been 
ascribed previously as the source of the significant colour gradients in spiral galaxies (|de Jong 
19961 : iGadotti & dos AnjoslEoOll : La Barbera et alJb003t iMacArthur et al.ll20o4. while ou r claim 



that the implied age gradients in spirals are small is only matched by iPeletier & Balcellsl (|1996|) . 
In the median, the metallicity gradients for all Virgo spiral types appear to be uniformly negative, 
while their age gradients are more varied. This complex pattern of Virgo spirals' colour gradients 
thwarts any realistic estimate of typical ages and metallicities for these galaxies as well. However , 
the colours of Virgo s pirals exhibit a morphological trend, also reported elsewhere (|de Jonglll996i 
Mac Arthur et al.l |2004|) . which may be explained by the earlier galaxies being older and more 
enriched than the later ones. 



3.2.2. Structure 

We show in Figure \\0\ the same colour-colour diagram as in Figure [7] (including the same 
model) but with the points differentiated according to the following structural parameters: C28, 
m H , [i e or r e . In these plots, we focus on galaxy centers (i.e., within an aperture of 0.5 r e ) to ease 
the distinction of colour trends as a function of galaxy structure. We also indicate with error bars 
the median values and rms dispersions of the colour distributions for all structural parameter bins. 

The Virgo galaxy central colours exhibit multiple trends against structural parameters. The 
trend versus C 2 g, suggests that the most concentrated Virgo galaxies are the oldest and most metal- 
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rich. Since C28 measures the shape/concentration of a galaxy's light profile, this trend essentially 
confirms that found above versus morphology between the Virgo spheroidal and disk galaxies. 
The other trends, with m H and fj, e , indicate that the more luminous and brighter Virgo galaxies 
tend to be older and more enriched than less luminous systems. Our detection of these trends 



ents reported in numerous studies 


(Vader et al. 


19881 


Peletier et al. 


1990; I 


Jalcells & Peletier! 


1994] 


de Jong 


996; Jansen et al. 


2000; 


Jerien et al. 


2000; 


Tamura & O - 


ita 


2003 


: Mac Arthur et al. 


2004; 


Liu et al 


2009) (but see Kim & Ann 


1990 


and La Barbera et al. 


2004 


for arguments against the 



luminosity dependence of gas-poor galaxies' colours). No significant colour trend is seen against 
r e other than the most compact Virgo galaxies being older and more metal-rich than larger ones. 
This nu ll result agrees with the w e ak size trends i n spiral galaxies' colour gra dients seen else- 
where (IMunoz-Mateos et all 120071 : Liu et all 120091) . though Larodi et all (120021) reported a disk 
scale length dependence in the colour gradients of field and group dlrr's. 



3.2.3. Environment 



Figure ITTlis the environmental corollary to Figure [TO] To investigate possible colour-environment 
trends, the Mil galaxies have been assigned in each of these figures to one of four equally- 
populated bins in the following parameters: D M&1 , S, or Def HI . We find from these plots that 
the median colours of the different D M&1 and £ bins all roughly coincide, suggesting that cluster- 
centric location of, or loc al density about, a Vi rgo galaxy has little influenc e on its stellar content. 
This result is matched by IKoo et al J 120051) and lMunoz-Mateos et al.1 (|2007l) who found no change 
in th e colours or colo ur gradients of gas-rich galaxies with local density. The contradictory claim 
by (IParodi et al. Il2002l) about a density dependence of galaxy colours is diffused by considering that 
those authors used a sample which spanned a wider density baseline (fields-cluster) than ours. 

Amongst Virgo gas-rich galaxies, we find a clean trend of colours with Defm, such that more 
deficient systems are older and more metal-rich than less deficient ones. However, the comparable 
r-H colours found amongst intermediate-De/ ff/ galaxies suggest that age effects dominate over 
those of metallicity in creating the observed trend. To our knowledge, this is the first report of 
a gas deficiency-colour correlation for Virgo galaxies. However, a latent correlation between the 
Defm measurements for our galaxies and their morphology, with a Pearson correlation coefficient 
of -0.68, makes it unclear whether the colour-.De/#/ trend is fundamental or not. 
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3.2.4. Dust 



It would also be prudent to address the effects of dust on the colour properties of Virgo cluster 
galaxies. The presence of dust in star-forming systems is of concern in stellar population studies 
since reddened broadband colours may be erroneously attributed to their stellar content. For in- 
stance, the colour gradients of a galaxy with a smoothly-varying dust distribution may be taken 
as due to combined negative age and metallicity gradients within it, a consequence of the three- 
fold age-metallicity -dust degeneracy. Even fac e-on systems may be suscep tible to such effects 
(IPeletier et al.lll999l M04; Paper II). Moreover. iGoudfrooij & de Jongl (|1995|) have shown that the 
shallo w optical colour gr adients of gas-poor giants may be explained by a diffuse dust component, 
while iGanda et al.1 (120091) have found evidence for severe (but localized) dust attenuation in the 
optical-NIR profiles of nearby spirals. 

In Figure [121 we assess the dust reddening of the central colours of Virgo spirals in the r-H 
versus g-r plane, where individual galaxies are identified according to their NIR axis ratio, b/a. In 
so doing, we find large dispersions in the colour distributions for edge-on systems, which suggests 
that some of them likely suffer from significant reddening. The colour centroids of the highest 
and lowest inclination spirals are nearly identical, however, which is inconsistent with dust being 
the main driver of these galaxies' central colours. That is, in a statistical sense, the effect of dust 
on our galaxy colours may be neglected. This conclusion finds ample support in the liter ature 
Jde Jong! 1 199ft IPeletier & Balcellslll996l: iGadotti & dos AnjosllioOll : iMac Arthur et~al]l2004l) . but 
does not preclude a possible contribution of dust to the colour gradients of individual systems (e.g., 
Peletier et al.lll995l : lLa Barbera et al.ll2003|) . Indeed, we find dramatic colour variations between 
adjacent radial bins in a few Virgo spirals, indicative of significant reddening at those (redward) 
locations (i.e., those colour variations would otherwise imply radical changes in local age and 
metallicity), but the incidence of such bins is small. 



4. CONCLUSIONS 

By combining high-quality optical (gri) and near-infrared (H) imaging for ~300 Virgo cluster 
galaxies, we have studied the resolved broadband colours of a representative galaxy sample span- 
ning the full range of structural parameters. Below is a summary of our main findings about Virgo 
galaxies' colour gradients and trends in their central colours with structural and environmental 
parameters: 

• In the median, the colour gradients of Virgo galaxies systematically decrease with morpho- 
logical type, such that the gradients of giant early-types are nearly flat, while those of later 
types are strongly negative (i.e., red centers to blue outskirts). The colour gradients of early- 
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type dwarfs are more negative than those for the giants; 

• The colour gradients of all Virgo galaxies principally arise from stellar metallicity gradients, 
while stellar age gradients contribute only to the colours of Virgo irregulars; 

• The central colours of Virgo galaxies correlate with concentration, luminosity, and surface 
brightness, such that galaxies become redder with an increase in each of these parameters; 

• While the central colours of Virgo galaxies do not depend on cluster-centric distance or 
local galaxy density, those of Virgo gas-rich galaxies correlate strongly with their neutral 
gas deficiency; 

• The central colours of Virgo spiral galaxies are uncorrelated with axial ratio, thus limiting 
the role of dust in setting their colours. 

While the above results reflect on the stellar content of Virgo galaxies and their structural and 
environmental dependencies, we can explicitly apply stellar population models to galaxy colours. 
This is what we do in Paper II as a way to constrain formation scenarios for the major galaxy types 
in the Virgo cluster. 

This research has made use of the NASA/IPAC extragalactic and GOLDMine databases. 
Funding for the SDSS has been provided by the Alfred P. Sloan Foundation, the Participating 
Institutions, the National Science Foundation, the US Department of Energy, the National Aero- 
nautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and 
the Higher Education Funding Council for England. The SDSS Web Site is http://www.sdss.org/ 
The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. 
The Participating Institutions are the American Museum of Natural History, Astrophysical Institute 
Potsdam, University of Basel, Cambridge University, Case Western Reserve University, University 
of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation 
Group, The Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Insti- 
tute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy 
of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astron- 
omy, the Max-Planck-Institute for Astrophysics, New Mexico State University, The Ohio State 
University, University of Pittsburgh, University of Portsmouth, Princeton University, the United 
States Naval Observatory, and the University of Washington. J.R. and S.C. acknowledge finan- 
cial support from the National Science and Engineering Council of Canada through a postgraduate 
scholarship and a Discovery Grant, respectively. We also thank Brent Tully for his help in the 
acquisition of //-band imaging of Virgo cluster galaxies (as reported in Ml 1). 



-17- 



REFERENCES 

Abraham, R. G., Ellis, R. S., Fabian, A. C, Tanvir, N. R., & Glazebrook, K. 1999, MNRAS, 303, 
641 

Adelman-McCarthy, J. K., et al. 2007, ApJS, 172, 634 
Adelman-McCarthy, J. K., et al. 2008, ApJS, 175, 297 
Balcells, M., & Peletier, R. F. 1994, AJ, 107, 135 
Barazza, F. D., Binggeli, B., & Prugniel, P. 2001, A&A, 373, 12 
Barazza, F. D., Binggeli, B., & Jerjen, H. 2003, A&A, 407, 121 
Bell, E. F, & de Jong, R. S. 2000, MNRAS, 312, 497 

Bergvall, N., Ronnback, J., Masegosa, J., &Oumlstlin, G. 1999, A&A, 341, 697 

Binggeli, B., Sandage, A., & Tammann, G. A. 1985, AJ, 90, 1681 

Boselli, A., Boissier, S., Cortese, L., & Gavazzi, G. 2008, ApJ, 674, 742 

Bothun, G. D., & Gregg, M. D. 1990, ApJ, 350, 73 

Bremnes, T., Binggeli, B., & Prugniel, P. 1998, A&AS, 129, 313 

Bremnes, T., Binggeli, B., & Prugniel, P. 1999, A&AS, 137, 337 

Bremnes, T., Binggeli, B., & Prugniel, P. 2000, A&AS, 141, 211 

Burstein, D., Faber, S. M, Gaskell, C. M, & Krumm, N. 1984, ApJ, 287, 586 

Cardiel, N., Gorgas, J., Sanchez-Blazquez, P., Cenarro, A. J., Pedraz, S., Bruzual, G., & Klement, 
J. 2003, A&A, 409,511 

Courteau, S., de Jong, R. S., & Broeils, A. H. 1996, ApJ, 457, L73 

de Jong, R. S. 1996, A&A, 313, 377 

Dressier, A. 1980, ApJ, 236, 351 

Faber, S. M. 1973, ApJ, 179, 731 

Ferreras, I., Lisker, T., Carollo, C. M., Lilly, S. J., & Mobasher, B. 2005, ApJ, 635, 243 
Ferreras, L, Lisker, T., Pasquali, A., & Kaviraj, S. 2009, MNRAS, 395, 554 



-18- 



Franx, M., Illingworth, G., & Heckman, T. 1989, AJ, 98, 538 
Gadotti, D. A., & dos Anjos, S. 2001, AJ, 122, 1298 
Gallagher, J. S., Ill, & Hunter, D. A. 1986, AJ, 92, 557 

Gallazzi, A., Chariot, S., Brinchmann, J., White, S. D. M., & Tremonti, C. A. 2005, MNRAS, 362, 
41 

Ganda, K., Peletier, R. R, Balcells, M., & Falcon-Barroso, J. 2009, MNRAS, 395, 1669 
Gavazzi, G, Bonfanti, C., Sanvito, G, Boselli, A., & Scodeggio, M. 2002, ApJ, 576, 135 
Gavazzi, G, Boselli, A., Donati, A., Franzetti, R, & Scodeggio, M. 2003, A&A, 400, 451 
Gavazzi, G, Boselli, A., van Driel, W., & O'Neil, K. 2005, A&A, 429, 439 
Gil de Paz, A., & Madore, B. F 2005, ApJS, 156, 345 

Gonzalez-Perez, V., Castander, F J., & Kauffmann, G 2011, MNRAS, 411, 1151 

Goudfrooij, P., Hansen, L., Jorgensen, H. E., Norgaard-Nielsen, H. U., de Jong, T., & van den 
Hoek, L. B. 1994, A&AS, 104, 179 

Goudfrooij, P., & de Jong, T. 1995, A&A, 298, 784 

Grant, N. I., Kuipers, J. A., & Phillipps, S. 2005, MNRAS, 363, 1019 

Hunter, D. A., & Elmegreen, B. G 2006, ApJS, 162, 49 

James, P. A. 1994, MNRAS, 269, 176 

James, P. A., Salaris, M., Davies, J. L, Phillipps, S., & Cassisi, S. 2006, MNRAS, 367, 339 

Jansen, R. A., Franx, M., Fabricant, D., & Caldwell, N. 2000, ApJS, 126, 271 

Jerjen, H., Binggeli, B., & Freeman, K. C. 2000, AJ, 119, 593 

Kim, K. O., & Ann, H. B. 1990, Journal of Korean Astronomical Society, 23, 43 

Ko, J., & Im, M. 2005, Journal of Korean Astronomical Society, 38, 149 

Kobayashi, C, & Arimoto, N. 1999, ApJ, 527, 573 

Koo, D. C, et al. 2005, ApJS, 157, 175 

Koopmann, R. A., & Kenney, J. D. P. 2004, ApJ, 613, 866 



- 19- 

La Barbera, R, Busarello, G., Merluzzi, P., Massarotti, M., & Capaccioli, M. 2002, ApJ, 571, 790 

La Barbera, R, Busarello, G., Massarotti, M., Merluzzi, P., & Mercurio, A. 2003, A&A, 409, 21 

La Barbera, R, Merluzzi, P., Busarello, G., Massarotti, M., & Mercurio, A. 2004, A&A, 425, 797 

La Barbera, R, de Carvalho, R. R., Gal, R. R., Busarello, G., Merluzzi, P., Capaccioli, M., & 
Djorgovski, S. G. 2005, ApJ, 626, L19 

La Barbera, R, & de Carvalho, R. R. 2009, ApJ, 699, L76 

La Barbera, R, de Carvalho, R. R., de la Rosa, I. G., Gal, R. R., Swindle, R., & Lopes, P. A. A. 
2010. larXiv: 1006.40361 

Liu, C.-Z., Shen, S.-Y, Shao, Z.-Y., Chang, R.-X., Hou, J.-L., Yin, J., & Yang, D.-W. 2009, 
Research in Astronomy and Astrophysics, 9, 1119 

MacArthur, L. A., Courteau, S., Bell, E., & Holtzman, J. A. 2004, ApJS, 152, 175 

MacArthur, L. A. 2005, ApJ, 623, 795 

MacArthur, L. A., Gonzalez, J. J., & Courteau, S. 2009, MNRAS, 395, 28 

McDonald, M., Courteau, S., & Tully, R. B. 2009, MNRAS, 394, 2022 

McDonald, M., Courteau, S., Tully, R. B., & Roediger, J. C. 2011, MNRAS, in press 

McLaughlin, D. E. 1999, ApJ, 512, L9 

Mei, S., et al. 2007, ApJ, 655, 144 

Menanteau, R, et al. 2004, ApJ, 612, 202 

Michard, R. 1999, A&AS, 137, 245 

Michard, R. 2000, A&A, 360, 85 

Moriondo, G., et al. 2001, A&A, 370, 881 

Moth, P., & Elston, R. J. 2002, AJ, 124, 1886 

Munoz-Mateos, J. C, Gil de Paz, A., Boissier, S., Zamorano, J., Jarrett, T., Gallego, J., & Madore, 
B. F. 2007, ApJ, 658, 1006 

Papaderos, P., Loose, H.-H., Thuan, T. X., & Fricke, K. J. 1996, A&AS, 120, 207 



-20- 

Parodi, B. R., Barazza, F. D., & Binggeli, B. 2002, A&A, 388, 29 
Patterson, R. J., & Thuan, T. X. 1996, ApJS, 107, 103 
Peletier, R. R, & Valentijn, E. A. 1989, Ap&SS, 156, 127 

Peletier, R. R, Davies, R. L., Illingworth, G. D., Davis, L. E., & Cawson, M. 1990, AJ, 100, 1091 

Peletier, R. R, Valentijn, E. A., Moorwood, A. R M., Freudling, W., Knapen, J. H., & Beckman, 
J. E. 1995, A&A, 300, LI 

Peletier, R. R, & Balcells, M. 1996, AJ, 1 11, 2238 

Peletier, R. R, Balcells, M., Davies, R. L., Andredakis, Y., Vazdekis, A., Burkert, A., & Prada, R 
1999, MNRAS, 310, 703 

Peletier, R. R, et al. 2007, MNRAS, 379, 445 

Pierini, D. 2002, MNRAS, 330, 997 

Searle, L., Sargent, W. L. W., & Bagnuolo, W. G. 1973, ApJ, 179, 427 
Silva, D. R., & Elston, R. 1994, ApJ, 428, 511 
Skrutskie, M. R, et al. 2006, AJ, 131, 1 163 

Solanes, J. M., Sanchis, T., Salvador-Sole, E., Giovanelli, R., & Haynes, M. P. 2002, AJ, 124, 2440 
Suh, H., Jeong, H., Oh, K., Yi, S. K., Ferreras, L, & Schawinski, K. 2010, ApJS, 187, 374 
Tamura, N, & Ohta, K. 2000, AJ, 120, 533 

Tamura, N, Kobayashi, C., Arimoto, N, Kodama, T., & Ohta, K. 2000, AJ, 1 19, 2134 
Tamura, N, & Ohta, K. 2003, AJ, 126, 596 

Taylor, V. A., Jansen, R. A., Windhorst, R. A., Odewahn, S. C., & Hibbard, J. E. 2005, ApJ, 630, 
784 

Terlevich, R., Davies, R. L., Faber, S. M., & Burstein, D. 1981, MNRAS, 196, 381 

Terndrup, D. M., Davies, R. L., Frogel, J. A., Depoy, D. L., & Wells, L. A. 1994, ApJ, 432, 518 

Tinsley, B. M., & Gunn, J. E. 1976, ApJ, 203, 52 



Tortora, C, Napolitano, N. R., Cardone, V. R, Capaccioli, M., Jetzer, P., & Molinaro, R. 2010, 
MNRAS, 407, 144 

Trager, S. C, Faber, S. M, Worthey, G., & Gonzalez, J. J. 2000, AJ, 1 19, 1645 
Tremonti, C. A., et al. 2004, ApJ, 613, 898 
Tamm, A., & Tenjes, P. 2006, A&A, 449, 67 

Tully, R. B., Verheijen, M. A. W., Pierce, M. J., Huang, J.-S., & Wainscoat, R. J. 1996, AJ, 112, 
2471 

Vader, J. P., Vigroux, L., Lachieze-Rey, M., & Souviron, J. 1988, A&A, 203, 217 
van Zee, L., Barton, E. J., & Skillman, E. D. 2004, AJ, 128, 2797 

West, A. A., Garcia-Appadoo, D. A., Dalcanton, J. J., Disney, M. J., Rockosi, C. M., & Ivezic, Z. 
2009, AJ, 138, 796 

Worthey, G., Faber, S. M, & Gonzalez, J. J. 1992, ApJ, 398, 69 

Worthey, G. 1994, ApJS, 95, 107 

Wu, H., Shao, Z., Mo, H. J., Xia, X., & Deng, Z. 2005, ApJ, 622, 244 
Zaritsky, D., Kennicutt, R. C., Jr., & Huchra, J. P. 1994, ApJ, 420, 87 



This preprint was prepared with the AAS L?Tr^C macros v5.2. 



-22- 



Table 1 . Median g-H colours and colour gradients of Virgo galaxies 



Morphology 




dig-H)* 
dr 




(g-H)o 


N 




r/r e (magr/) 


r (mag kpcr') 


log r (mag dex 1 ) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


dSO 


-0.126 ±0.199 


-0.080 ±0.165 


-0.198 ±0.290 


3.000 ± 0.220 


19 


dE 


-0.177 ±0.208 


-0.167 ±0.196 


-0.309 ± 0.377 


2.843 ± 0.184 


51 


E 


-0.084 ±0.199 


-0.075 ± 0.339 


-0.267 ± 0.572 


3.253 ± 0.353 


32 


SO 


-0.046 ± 0.212 


-0.028 ±0.152 


-0.115 ±0.310 


3.283 ± 0.280 


53 


Sa-Sb 


-0.138 ±0.314 


-0.057 ± 0.080 


-0.322 ± 0.424 


3.274 ± 0.578 


24 


Sbc+Sc 


-0.190 ±0.324 


-0.066 ±0.101 


-0.428 ± 0.407 


3.176 ±0.663 


15 


Scd+Sd 


-0.259 ± 0.412 


-0.106 ±0.180 


-0.371 ±0.652 


3.037 ± 0.588 


17 


Sdm+Sm 


-0.017 ± 0.467 


-0.009 ±0.231 


-0.072 ± 0.616 


2.459 ± 0.420 


9 


Im 


-0.181 ±0.391 


-0.076 ± 0.278 


-0.228 ± 0.662 


2.302 ± 0.499 


13 


BCD 


+0.051 ±0.233 


+0.071 ± 0.321 


+0.197 ±0.403 


2.355 ± 0.417 


8 


s? 


+0.141 ±0.400 


+0.108 ±0.213 


+0.158 ±0.468 


2.415 ± 0.466 


9 


? 


-0.137 ±0.178 


-0.077 ± 0.266 


-0.246 ± 0.479 


2.850 ± 0.422 


8 



"The median colour gradients for Virgo galaxies were computed for three choices of radial coordinate: 
effective radii (left), kiloparsec (center), and logarithmic (right). 
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Fig. 1. — (left) Comparison of Defm measurements for Virgo galaxies in common between S02 
and KK04. The red line has slope unity, while the solid and dashed black lines represent a linear 
fit to the data and the la deviation about the fit, respectively. The Pearson (linear) correlation 
coefficient r and a of the fit are both provided in the top-left corner, (right) Same as left, but 
comparing S02 and G05. 
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Fig. 2. — g-// profiles as a function of scaled radius, r/r e , for Virgo galaxies binned by morpho- 
logical type. The median profiles and their la-dispersion are indicated by the solid and dashed 
lines, respectively. 
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Fig. 3. — Size of the first radial bin versus the FWHM of the maximum seeing disk for Virgo 
galaxies. The 1 : 1 relation is marked by the red line. 
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Fig. 4. — Comparison of the g-H colour gradients for Virgo galaxies determined using either a 
fixed or variable radial range in the linear fits to the colour profiles of individual galaxies. The data 
point colours reflect the maximum extent of each galaxy colour profile. For the fixed radial range 
fits, only data points within r < 2r e were included. The black line marks the 1 : 1 relation. 
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Fig. 5. — Comparison of the median g-H colour gradient estimates (as determined in r e space) for 
Virgo galaxies, binned by morphological type, using a minimum S/N of either 10 or 50. The size 
of the error bars reflects the 1 a (rms) dispersion within each morphological bin. 
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Fig. 6. — The (potential) contribution of sky subtraction errors to the estimated colour gradients 
(in r e space) of individual Virgo galaxies, binned by morphological type. The red line connects the 
median error within each morphological bin. 
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Fig. 7. — (left) r-H versus g-r colour-colour diagrams for Virgo dE and dSO systems. Stars and 
filled circles represent the radial bins which enclose the galaxy centers and effective radii, respec- 
tively, while the coloured lines show the median colour profiles for both galaxy types. The grid 
represents a stellar population model for an exponential star formation history with no chemical 
evolution. The dashed and solid black lines correspond to lines of constant mean age and metal- 
licity, respectively, with representative values indicated throughout the grid. The metallicities are 
expressed in percentages, according to mass contribution (i.e., Z = 0.02). The error bars shown 
at lower-right indicate the typical uncertainties in the measured colours due to both calibration 
and sky subtraction errors for the entire sample (the latter is shown for both galaxies' central and 
outermost radial bins). The orange-red arrow in the top-left corner shows the reddening vector of 
a foreground dust screen model for the indicated amount of extinction, (right) Same as left but for 
Virgo E galaxies. 
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Fig. 8. — Same as Fig. [7]for Virgo lenticular {left) and irregular {right) galaxies. 
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Fig. 10. — (top left) Variation in the central (r < 0.5r e ) r-H and g-r colours with concentration 
C28 for Virgo galaxies. Galaxies have been grouped into equally-populated bins by C28, with the 
bin ranges and point colours indicated in the lower-right corner of the window. The overplotted 
error bars show both the centroid and rms scatter of the colour distribution within each C28 bin. 
The other plots are the same as at top-left, but show the variations in central r-H and g-r colours 
of Virgo galaxies with apparent magnitude m H (top right), effective surface brightness fi e (bottom 
left), and effective radius r e (bottom right). 
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Fig. 1 1 . — Same as Fig. [[0]but showing variations in central r-H and g-r colours of Virgo galaxies 
with distance from M87 D M si (top left), local galaxian surface density E (top right) and neutral gas 
deficiency Def HI (bottom). The plot against Def HI only includes 101 galaxies due to the limited 
Hi mapping of the Virgo cluster. 
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Fig. 12. — Same as Fig. [I0]but showing variations in central r-H and g-r colours of Virgo galaxies 
with //-band axis ratio, b/a. Only the colours of Virgo spiral galaxies are plotted here. 



